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Abstract 

We observe a signal for the doubly charmed baryon H^ in the decay mode H^ -^ pD~^K~ 
to complement the previous reported decay H+ -^ A+ K~ir~^ in data from SELEX, 
the charm hadro-production experiment at Fermilab. In this new decay mode we 
observe an excess of 5.62 events over a combinatoric background estimated by event 
mixing to be 1.38 it 0.13 events. The mixed background has Gaussian statistics, 
giving a signal significance of 4.8cj. The Poisson probability that a background fluc- 
tuation can produce the apparent signal is less than 6.4 x 10~^. The observed mass 
of this state is 3518 it 3 MeV/c^, consistent with the published result. Averaging the 
two results gives a mass of 3518.7 it 1.7 MeV/c^. The observation of this new weak 
decay mode confirms the previous SELEX suggestion that this state is a double 
charm baryon. The relative branching ratio for these two modes is 0.36 it 0.21. 

Key words: Doubly Charmed Baryon 
PACS: 14.20.Lq, 14.40.Lb, 13.30.Eg 



1 Introduction 



In 2002, the SELEX collaboration reported the first observation of a candidate 
for a double charmed baryon, decaying as S+ — > A+ K~7r+[lj]. The state had 
a mass of 3519 ± 2 MeV/c^, and its observed width was consistent with exper- 
imental resolution, less than 5 MeV/c^. The final state contained a charmed 
hadron, a baryon, and negative strangeness (A+ and K~), consistent with the 
Cabibbo-allowed decay of a S^ configuration. In order to confirm the inter- 
pretation of this state as a double charm baryon, it is essential to observe the 
same state in some other way. Other experiments with large charm baryon 
samples, e.g., the FOCUS and E-791 fixed target charm experiments at Fer- 
milab or the B-factories, have not confirmed the double charm signal. This 
is consistent with the SELEX results. The report in Ref. [1| emphasized that 
this new state was produced by the baryon beams (S~, proton) in SELEX, 
but not by the vr" beam. It also noted that the apparent lifetime of the state 
was significantly shorter than that of the A^, which was not expected in a 
model calculation based on Heavy Quark Effective Theory J2|. 

Another way to confirm the S+ is to observe it in a different decay mode 
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that also involves a final state with baryon number and charm (not anti- 
charm). One such mode, involving only stable charged particles, is the channel 
H+ -^ pD^K~ . Observing a mass peak near the S+ -^ A+ K~ti^ peak at 
3519 MeV/c^reported in Ref. jl| in a channel combining a proton with a D"'"K~ 
pair but not a D^K^ pair would confirm the existence of the H^ state. Here 
we report the first observation of S^ -^ pD^K~. 



2 Experimental apparatus 



The SELEX experiment used the Fermilab charged hyperon beam at 600 GeV 
to produce charm particles in a set of thin foil targets of Cu or diamond. 
The negative beam composition was about 50% S~, 50% vr". The positive 
beam was 90% protons. A beam Transition Radiation Detector identified each 
beam particle as meson or baryon with zero overlap. The three-stage magnetic 
spectrometer is shown elsewhere J3|, y]. The most important features are the 
high-precision, highly redundant, vertex detector that provided an average 
proper time resolution of 20 fs for the charm decays, a 10 m long Ring- Imaging 
Cerenkov (RICH) detector that separated vr from K up to 165 GeV/c [5], and a 
high-resolution tracking system that had momentum resolution of (Jp/p < 1% 
for a 150 GeV/c proton. 

The experiment selected charm candidate events using an online secondary 
vertex algorithm. A scintillator trigger demanded an inelastic collision with at 
least four charged tracks in the interaction scintillators and at least two hits 
in the positive particle hodoscope after the second analyzing magnet. Event 
selection in the online filter required full track reconstruction for measured fast 
tracks (j9>15 GeV/c). These tracks were extrapolated back into the vertex 
silicon planes and linked to silicon hits. The beam track was measured in 
upstream silicon detectors. A three-dimensional vertex fit was then performed. 
An event was written to tape if all the fast tracks in the event were inconsistent 
with having come from a single primary vertex. This filter passed 1/8 of all 
interaction triggers and had about 50% efficiency for otherwise accepted charm 
decays. The experiment recorded data from 15.2 x 10^ inelastic interactions 
and wrote 1 x 10^ events to tape using both positive and negative beams. The 
sample was 65% S~-induced, with the balance split roughly equally between 
7r~ and protons. 

The offline analysis selected single charm events with a topological identifica- 
tion procedure. Only charged tracks with reconstructed momenta were used. 
Tracks which traversed the RICH {p > 22 GeV/c) were identified as protons or 
kaons if those hypotheses were more likely than the pion hypothesis. All other 
tracks were assumed to be pions. The primary vertex was refit offline using 
all found tracks. An event was rejected if all tracks were consistent with one 



primary vertex. For those events which were inconsistent with a single primary 
vertex, secondary vertices were formed geometrically and then tested against 
a set of charge, RICH identification, and mass conditions to identify candi- 
dates for the different single charm states. Candidate events were written to 
a charm data summary file. Subsequent analysis began by selecting particular 
single-charm species from that set of events. 



3 Search Strategy 

In this study we began with the SELEX D^ sample that has been used in life- 
time and hadropro duct ion studies J7| . The sample-defining cuts are described 
in that reference. No new cuts on the D mesons were introduced in this anal- 
ysis. The D meson momentum vector had to point back to the primary vertex 
with x^ < 12 (the double charm lifetime is known to be much shorter than 
the D meson lifetime, so the D meson pointback is not affected by having come 
from a secondary decay). The D meson decay point must have a vertex separa- 
tion significance of at least 10 a from the primary. Everywhere in these analyses 
the vertex error used is the quadrature sum of the errors on the primary and 
secondary vertices. The K was positively identified by the RICH detector. The 
pions were required to be RICH-identified if they went into its acceptance. The 
D^ — >■ K'TT^TT^ and D'^ -^ K^t^~tt~ mass distributions are shown in Fig. 1. 
There are 1450 D^ -^ K^ti^tt^ decays and 2450 D~ -^ K^tt'tt' decays in 
these samples. The D^ -^ K^-k^-k^ events contribute to the signal channel. 
The D^ -^ K^7i^7T~ events cannot come from the decay of a double charm 
baryon and will be used as a topological background control sample. The yield 
asymmetry stems from the d-quark contribution of the S~ beam component 
that gives a sizeable production asymmetry favoring leading D~ production 
over D^ production, as we have reported for other charm systems [y]. 

The track-based search code is identical to that used on the A+ sample in the 
original investigation [l[ . The premise is that a ccd state will make a secondary 
decay vertex between the primary production vertex in one of the thin foil tar- 
gets and the observed D meson decay vertex, which must lie outside material. 
We looked for intermediate vertices using all charge zero pairs of tracks from 
the set of reconstructed tracks not assigned to the D-meson candidate. The 
additional positive track in this final state must be RICH-identified as a pro- 
ton if it traverses the RICH. The negative track in the new vertex is assigned 
the kaon mass. (This track typically missed the RICH acceptance by being 
too soft or too wide-angle, as confirmed by simulation.) We made background 
studies by (i) assigning the negative track a pion mass, (ii) looking for proton 
plus positive track combinations with a D meson, and (iii) looking at proton 
D~K~ combinations (wrong-sign charm). We require a good 3-prong vertex 
fit with a separation significance of at least 1.0 cr from the primary vertex. 
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Fig. 1. D^ —* K TT^TT^ (left) and D -^ K^tt it (right) mass distributions with 
cuts used in this analysis 

the same requirement used in Ref . [l| . The primary position was recalculated 
from the beam track and secondary tracks assigned to neither the D nor the 
pK~ vertices. Results presented in this paper come from this analysis. 



4 Search results and significance 



The signal search was based on a 10 MeV/c^window centered on the S+ 
mass 3519 MeV/c^ from Ref. [l|. The expected mass resolution for the de- 
cay S+ -^ pD^K~ is 4 MeV/c^. Our simulation correctly reproduces the 
observed widths of all our reported single charm mesons and baryons. The 10 
MeV/c^ window should collect 80% of the events in a S+ signal. The results 
are insensitive to changing bin boundaries by up to half a bin. The background, 
assumed to be fiat, is evaluated outside a 20 MeV/c^ window, to avoid putting 
the remaining 20% of any signal into the background. 

The right-sign mass combinations in Fig. 2(a) show an excess of 5.4 events 
over a background of 1.6 events. The wrong-sign mass combinations (c quark 
in the decay) for the pD^n^ final state are also plotted in Fig. 2 (b), scaled 
by 0.6 for the D+/D^ratio. The wrong-sign background shows no evidence 
for a significant narrow structure near 3519 MeV/c^. The average wrong-sign 
occupancy is 0.4 events/bin, exactly the background seen in the right-sign 
channel. This confirms the combinatoric character of the background popula- 
tion in the right-sign signal. We have investigated all possible permutations of 
particle assignments. The only significant structure observed is in the channel 
S+ —>■ pD^K~ , the place where a double charm baryon decay can occur. 

In order to assign a significance level to this peak, we have combined statisti- 



4 



^ 3.5 

c 

LU 3 



(a) pD^K" 



peak mass: 
3516 MeV 



4-bin Poisson Prob 
< 6.4 X 10"^ 
L/o> 1.0 



2.5 
2 

1.5 

1 hn n 
0.5 



•k ' 

I 1-75 



1.5 



1.25 



c 

s 

2 0.75 
CO 

0.5 



0.25 



146 3.48 



3.5 3.52 3.54 



3.56 3.58 
M(p D* K) 



(b) wrong sign: p D" K" 
L/a>1.0 



3.46 3.48 3.5 3.52 3.54 3.56 3.58 

M(p D K) 



Fig. 2. (a) H^ -^ pD^K~ mass distribution for right-sign mass combinations. Ver- 
tical dashed lines indicate the region of smallest fluctuation probability as described 
in the text, (b) wrong-sign events with a pD~K"'", scaled by 0.6 as described in the 
text. The horizontal line shows a maximum likelihood fit to the occupancy. 



cal methods used in the original double charm paper [1] and the event-mixing 
method used in Ref. J8|. We set out to test the hypothesis that the back- 
ground events in Fig. 2 are random combinatoric tracks associated with real 
D+ mesons. To mix events we took a D"*" meson in the peak region (Fig. 1) and 
combined it with proton and K~ tracks extracted from other events. Each D"*" 
was reused 25 times. To compare to the combinatoric background in Fig. 2, 
we scale the mixed-event background down for the multiple D+ usage. 



The resulting background distribution predicts the observed distribution very 
well. The mean number of background events below, in, and above the 4-bin 
signal peak is 5.74 ± 0.26, 1.38 ± 0.13, and 10.60 ± 0.36. This agrees well with 
the 8, 1.6, and 10 events that we observe in the corresponding regions. A mass 
plot with the combinatoric background level is shown in Fig. 3. The back- 
ground that we observe is completely consistent in shape and normalization 
with random combinatoric tracks associated with real D"*" mesons. 
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Fig. 3. 



pD^K mass distribution from Fig. 2(a) with high-statistics mea- 



surement of random combinatoric background computed from event-mixing. 

We take the combinatoric background to give a proper measure of the expected 
background under the 7-event signal region. The single-bin significance of the 
signal, using the above background numbers, is: S = (7 - 1.38)/-\/l-38 -|- 0.13^ 
= 4.8 a. The Poisson probability of observing at least this much excess, includ- 
ing the Gaussian uncertainty in the background, is 6.4 x 10"'^. Both of these 
statistical significance calculations use methods identical to those in Ref [l|. 
This indicates a robust signal atop a combinatoric background whose shape 
and normalization are very well understood. 



5 Signal Properties 



In order to estimate the mass of the H^ -^ pD~^K~ state in light of the sparse 
statistics in Fig. 2, we fixed the width of the Gaussian to 4 MeV/c^and fitted 
the data distribution around the signal peak. The fit mass is 3518± 3MeV/c^. 



This agrees beautifully with the measurement of 3519 ± 2 MeV/c^from the 
original double charm baryon report. We present these data as confirmation of 
the double charm state at 3520 MeV/c^ in a new decay mode S+ — >■ pD~^K~ . 
The weighted average mass is 3518.7 ± 1.7 MeV/c^. The mass distributions 
for the two channels are shown in Fig. 4. 
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We have used the simulation to study the relative acceptance for the two de- 
cay channels S+ — >• pD^K~ and S+ — >• A+ K~7i^ in order to quote a relative 
branching ratio. The overall acceptance, including the single charm selection 
and the proton ID requirements in the S+ -^ pD~^K~ mode, is very similar. 
SELEX measures the relative branching ratio F(S+ — > pD^K^)/T{E'^^ — > A+ 
= 0.36±0.21. The systematic error due to acceptances is well understood from 
single charm studies and is negligible compared to the statistical error. 

In Ref. pj we noted that all observed ccd events were produced by the baryon 
beams. None came from pions. In this sample, 1 event out of the 7 in the peak 
region seen in Fig. 2 is a pion beam event, and 1 of the 19 sideband events 
comes from the pion beam sample. This sample is consistent with the view 
that double charm baryons are produced dominantly by the baryon beams 
in SELEX. In another comparison, we had noted that the S^ -^ A+ K~7t^ 
decays had an exceptionally short reduced proper time distribution, indicating 
a S+ decay lifetime 5-10 times shorter than the A+ lifetime. That feature is 
confirmed by the S+ -^ pD^K~ channel. As we noted in Ref. [l|, our lifetime 
resolution is excellent but we cannot exclude lifetime (strong decay) for these 
events. The width of this peak is completely consistent with simulation of a 
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zero-width state, unlikely for a strong decay of a massive state. Also, we do 
not see an increase in the signal when we reduce the vertex significance cut 
L/(j below 1. If this were a strong decay, one would expect as many events 
with L/cr of -1 as +1, so the signal should grow significantly. It does not. 

In Ref. yj we noted that the S+ — ^ A+ K'-n^ yield and acceptance implied 
that a large fraction of the A+ decays seen in SELEX came from double charm 
decays. That was a surprise. For the S+ — > pD~^K~ case that is not true. Only 
a few percent of the SELEX D"*" events are associated with double charm. 



6 Summary 



In summary, SELEX reports an independent confirmation of the double charm 
baryon H+ previously seen in the H^ — > A+ K~7t~^ decay mode via the obser- 
vation of its decay into the H+ -^ pD^K^ final state. Using only very loose 
cuts gives the statistically-significant signal shown in Fig. 2. A combinatoric 
background from event mixing describes the non-signal distribution well. The 
Gaussian signficance using this background estimate is 4.8a. This decay mode 
confirms that the S+ lifetime is very short and that it is produced dominantly 
by baryon beams, as we reported in Ref. [l|. 
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